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Abstract 

We study relativistic S+Au collisions at 200^4 GeV/c using a hydrodynam- 
ical approach. We test various equations of state (EOSs), which are used to 
describe the strongly interacting matter at densities attainable in the CERN- 
SPS heavy ion experiments. For each EOS, suitable initial conditions can 
be determined to reproduce the experimental hadron spectra; this empha- 
sizes the ambiguity between the initial conditions and the EOS in such an 
approach. Simultaneously, we calculate the resulting thermal photon and di- 
electron spectra, and compare with experiments. If one allows the excitation 
of resonance states with increasing temperature, the electro-magnetic signals 
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from scenarios with and without phase transition are very similar and are not 
resolvable within the current experimental resolution. Only EOSs with a few 
degrees of freedom up to very high temperatures can be ruled out presently. 
We deduce an upper bound of about 250 MeV for the initial temperature from 
the single photon spectra of WA80. With regard to the CERES dilepton data, 
none of the EOSs considered, in conjunction with the standard leading order 
dilepton rates, succeed in reproducing the observed excess of dileptons below 
the p peak. Our work, however, suggests that an improved measurement of 
the photon and dilepton spectra has the potential to strongly constrain the 
EOS. 

PACS numbers: 25.70.np, 12.38.Mh, 12.40.Ee, 47.75.+f 
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I. INTRODUCTION 



The use of hydrodynamics for simulations of nuclear collisions has a long tradition |T]-0 . 
If applicable, hydrodynamics has some advantages over the more fundamental kinetic cal- 
culations, which are usually performed as Monte-Carlo simulations. Besides its relative 
simplicity, the use of familiar concepts such as temperature, flow velocity, energy density, 
pressure, etc. leads to an intuitively transparent picture of the space-time evolution of the 
system. Another great advantage is the direct use of the equation of state (EOS) of strongly 
interacting matter. Testing different EOSs by comparing with experiments should give us 
more insight into the behavior of nuclear matter under extreme conditions of temperature 
and/or density. This is important, since one of the main reasons to study high energy 
heavy ion collisions is to confirm the possible phase transition from hadronic matter to the 
quark-gluon plasma (QGP) The necessary energy densities are estimated to be around 
1-2 GeV/im^, presently experimentally available at the Brookhaven AGS and the CERN- 
SPS. Of these two facilities, conditions to observe signals of the phase transition are more 
favorable at the CERN-SPS, because of the higher incident momenta of the nuclei. 

In this work, we analyze the data from the heavy ion experiments at the CERN-SPS 
with 200y4 GeV/c incident momentum [Q- concentrating on the S + Au system. For 
this system, yields of both the hadronic and electromagnetic (dileptons and photons) probes 
are now available. A hydrodynamical treatment of the nuclear collisions in this energy range 
is not new [|1^]- [|19[, but the large set of new and updated experimental data allows us to 



achieve a deeper insight into the space-time evolution of these reactions. 

Our emphasis here is on the simultaneous description of the hadron and electromagnetic 
data. Since hadrons interact throughout the dense phase of the collision, the hadronic 
spectra are chiefly determined from the conditions at the freeze-out of matter. However, 
different initial conditions combined with different EOSs can lead to the same final particle 
distributions. In contrast, photon and lepton pairs are emitted throughout the dense stage 



and escape without interactions. The measured spectra thus probe the different temperature 
and ffow conditions during the evolution. These should have a different dependence on the 
initial conditions and on the EOS, in comparison to the hadron spectra. This may help to 
further reduce the ambiguity in the initial conditions and uncertainties in the EOS. 

Our hydrodynamical calculations assume azimuthal symmetry around the beam axis. 
Thus, we simulate only central collisions with impact parameters close to zero. We do not 
describe the production of matter within hydrodynamics, but start the calculation at an 
initial time when the system is likely to have reached thermal and chemical equilibrium. 
We explore different initial conditions, taking guidance from the experimental data (where 
available) on particle production in proton-proton (p + p) collisions in choosing what could 
be considered as a realistic initial state. However, one must be aware that nuclear effects are 
likely already presentat the production stage, and the initial conditions cannot be uniquely 
determined from p + p processes alone. 

A major concern is the uncertainty in the nuclear EOS. It is experimentally and theoret- 
ically well known only around the ground state of nuclear matter. At higher temperatures 
and densities, model predictions are widely different. The main goal in this study is to derive 
constraints on the EOS by comparing the hydrodynamical calculations with hadronic and 
electromagnetic data simultaneously from S -|- Au coUisions at 200A GeV/c. 

The paper is organized as follows: In Sec II, we briefly introduce the hydrodynamical 
equations and address some problems which arise when seeking numerical solutions of these 
equations. In Sec. Ill, the different EOSs explored are described. The initial conditions are 
discussed in Sec. IV. In Sees. V, VI and VII, we describe the calculation of the final hadron, 
photon, and dilepton spectra, respectively. Results are discussed in Sec. VIII, and our 
conclusions are given in Sec. IX. In the appendix, we briefly explain how we incorporated 
the CERES kinematic cuts and the detector resolution. 
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II. THE HYDRODYNAMICAL EQUATIONS 



The basic equations of hydrodynamics are the local conservation of energy and momen- 
tum, which in Lorentz-covariant form are written as 

9^T^^(x) = . (1) 

We use the ideal fluid ansatz for the energy-momentum tensor 

r^-(a;) = [6{x) -p{x)g^'', (2) 

where e{x) is the local energy density, p{x) is the local pressure, and u^{x) is the local four- 
velocity, normalized to u^u^ = 1. In principle, viscous effects may be included, but would 



lead to a major increase in the numerical effort 0j2^. However, the calculation generates 
some numerical viscosity, as explained below. 

We include finite baryon density, ps, in our system and express baryon number conser- 
vation locally in the form 

d,J^ = , (3) 

in terms of the baryon current = PbU^- 

Using the definition Eq. (^, the continuity Eqs. (j^) and (0) can be written more explicitly 

as 

dtT°° + V ■ (T°V) = -V ■ (pv) , 

9iT0* + V-(T°V) = -9ip, 

dtfB + V ■ (jM = . (4) 

In order to solve Eq. (1), one additional input is needed, namely, the equation of state 
(EOS). The EOS relates the pressure p = p{e,pB) to the energy and baryon densities. The 
different choices of the EOS are described in the next section. 



The partial (hyperbolic) differential equations, Eq. (^, are solved numerically on a com- 
puter using a finite difference method. We use the SHASTA algorithm in two spatial 
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dimensions, radial and longitudinal. Some details of implementing this algorithm are ex- 



plained in the appendix of Ref. [|20[. However, there are two essential modifications. Firstly, 
since we now have a problem with two spatial dimensions |2^, the flux correction must 
be modified accordingly, as described in Ref. Secondly, the flux correction involves 

a parameter 77 [^,^, the anti-diffusive constant, which regulates the residual numerical 
anti-diffusive flux in the algorithm. Theoretically, for rj = 1/8, diffusion vanishes almost 
completely, leading, however, to the appearance of ripples in the presence of steep gradients. 

On the other hand, a residual diffusive flux mimics viscosity in the calculation, since it 
generates entropy in the system. In Ref. a large residual diffusive flux was allowed by 
choosing a small value for the anti-diffusive constant, 77 = 0.08. This value was determined 
by trying to create as much entropy as is maximally allowed in the rarefaction wave, typically 
< 10%. However, it turns out that the flux across the freeze-out surface (Eq. (|27|) below) now 
has two contributions; the normal convective part and the diffusive part. Because we cannot 
assign a velocity to the numerical diffusive flux, it cannot be converted into particle spectra 
at the freeze-out surface. Therefore, we diminish the diffusive flux by setting rj = 0.11, as 



suggested also in Ref. |2J|, and neglect the small residual diffusive flux across the freeze- 
out hypersurface. This leads to a small, of order <5%, loss in the baryon number and 
total energy across the freeze-out hypersurface. (The stability of the numerical calculation 
increases as rj is diminished from the limiting value 1/8 = 0.125. Our choice, r] = 0.11, is a 
compromise between numerical stability and the diffusive loss across the freeze-out surface.) 
The value of rj has a small influence on the space-time evolution, a larger diffusive flux 
leading to a more rapid cooling of the system. This has only a small effect, especially on 
the electromagnetic spectra. 

At the moment, our numerical code assumes a mirror symmetry in the longitudinal 
direction with respect to the center of mass {f{z) = f{—z) for scalar quantities). This is a 
very good approximation for symmetric colliding systems. However, we apply it here to the 
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asymmetric S + Au collisions. Thus, we are presently unable to reproduce the asymmetry 
in the rapidity distributions of hadrons, which is observed in experiment. However, in 
the central rapidity region, where most of the experimental data are measured, we do not 
expect this to affect significantly the transverse spectra or the electromagnetic yields. We 
will consider the longitudinal asymmetry in a later work. 



III. THE EQUATION OF STATE 

The equation of state (EOS) in the energy density domain of e>l GeV/fm^ or baryon 
density pe ^ 2po, where po is the density of ground state nuclear matter, is theoretically quite 
uncertain. At the CERN-SPS energies, the relevant energy densities are not high enough 
that the well established techniques of finite temperature perturbative QCD will apply. At 
the same time, the energy densities are typically too high for low energy approaches, such as 
chiral perturbation theory, to be applicable. Hadronic and quark-gluon matter at these tem- 
peratures and densities, and especially the chiral and/or de-confinement transitions between 
the two phases at some critical temperature and density, are non-perturbative phenomena 
which are currently investigated through lattice simulations of QCD. 

Simulations of pure SU(3) gauge theory indicate the occurence of a first-order deconfine- 
ment transition around ~ 260 MeV through studies of hysterisis, co-existing states, and 
abrupt quantitative changes in the various thermodynamic functions. Recent simulations 
on larger lattices (16^ x 32) suggest that finite size effects are reasonably well understood 
for pure glue — an estimate for the critical energy density for pure SU(3) is ~ 1 GeV/fm^ 



25| . For full lattice QCD with dynamical quarks, the situation is less clear, particularly in 



the critical region. However, a recent calculation |^ shows that the energy density displays 



a rapid change in a narrow region around Tc ~ 150 MeV, while the pressure changes more 
smoothly. 

Though the present uncertainities preclude a quantitative description of the thermody- 
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namic functions around the critical region, one may observe i) a smooth crossover between 
hadronic matter and the quark gluon phase and ii) a sizable and rapid change in the entropy 
density. These features may be roughly reproduced in a two phase description of the transi- 
tion. Indeed, even if the transition were a sharp first order transition, it is unlikely that the 
hydrodynamic fiow simulations would be extremely sensitive to the width (in temperature) 
of the critical region 1^^. Therefore, for the purpose of our hydrodynamic simulations, we 



will construct EOSs for both the hadron and quark gluon state variables and match them 
at the critical boundary in temperature and chemical potential by a Maxwell construction. 
In practice, when implemented in simulations, this boundary is smoothed over using a hy- 
perbolic tangent profile of width A. It has been checked that the results of our simulations 



are insensitive to A, as long as A/T^ << 1 [|28|. In addition to studying EOSs with a sharp 



cross over (but differing critical temperatures and latent heat), we will also consider an EOS 
with purely hadronic degrees of freedom. 

Let us first discuss the EOS we use for the hadronic phase. For a dilute hadronic gas 
at temperatures well below the pion mass, state variables can be computed reliably using a 
virial expansion with input from empirical scattering cross sections [p9| , p0| . It was observed 
that the state variables for an interacting hadron gas are well approximated by those for 
a Boltzmann gas of free hadrons and resonances. At temperatures comparable to the pion 
mass or higher, third and higher virial coefficients are important. Presently, there is no 
reliable way to compute these. For want of a systematic prescription, we shall stretch 
our conclusions from the virial expansion approach and assume that a dense hadron gas is 
roughly approximated by a gas of free hadrons and resonances. 

We restrict our studies here to a hadron gas, which contains the following hadrons with 
their corresponding anti-particles: 

71, K, r], p, u, K*, p, n, r]', (/), A, S, A, ai, S, S(1385). (5) 

Beside stable particles (on strong interaction time scales), the hadron gas also contains 



resonance states. This should mimic the attractive interaction among the hadrons in the 



spirit of the bootstrap model of Hagedorn 31 1, although we keep the number of resonance 



states finite. If we assume a phase transition to the QGP at temperatures around Tc ~ 150 
MeV, then the hmited number of hadronic states is justified, because of the suppression of 
higher mass states through the Boltzmann factor. 

As mentioned above, we will also consider a pure hadronic EOS, which does not exhibit 
a phase transition. In this case, somewhat higher temperatures are encountered in our 
calculations. Therefore, a sensitivity on the upper mass cut in the hadronic mass spectrum 
could be expected. We expect this to show up more in the values of the temperature rather 
than in the evolution of the flow. While hadronic spectra are not very sensitive to this 
cut, electromagnetic rates, however, could decrease. On the other hand, for consistency, 
we should include the increase in electromagnetic emission from processes involving the 
heavy resonances. To some extent, these effects tend to cancel each other. Mainly due to 
the numerical limitations of our hydrodynamical simulations, we have included only those 
states listed above, even at high temperatures. 

In order to derive the EOS, we begin with the grand canonical partition function for a 
non-interacting resonance gas. It is given by 

Z«(r, V, /iB, /is) = n exp [ZniT, V, fx,)] , (6) 

h 

where V is the volume, and the product is over the different hadron species h. The chemical 
potential fi, of the hadron h is given by its baryon number and its strangeness through 

Hh = BhUB + Shjxs (7) 

in full chemical equilibrium. We will use this assumption throughout the expansion, although 
there are indications of deviations, especially for the strange particles [^. The partition 
function Eq. is, in general, a function of four unknown variables. For calculations of 
intensive variables, the volume cancels and the value of the strangeness chemical potential 
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/is can be related to T and /ib by the requirement of local strangeness neutrality. 
The partition function Z^^T, V, fj,h) for hadron species h is 

In Z,,{T, V, /i,) = [3VpH = J dE ^^(^_^^) ^ ^ , (8) 

where Qh is the degeneracy factor, rrih is the mass, and /i^ is the chemical potential of hadron 
h. The ± sign corresponds to fermions and bosons, respectively. From the partition function 
Eq. (^, we can calculate all thermodynamical quantities. Specifically, we have 

1 dliiZ^ , , 

^(T.f-B) = ^— (10) 

p(r,^B) = r^ = f in2«. (11) 

From these equations, the pressure can be solved as a function of the energy density e 
and the baryon density pe- This form of the EOS, p = p(£:,pb), is needed in solving the 
hydrodynamic equations and in practice is obtained numerically from Eqs. (|T0|) and (]Tl|). 

In the limit of high baryon densities, the repulsive interactions between particles have 
to be taken into account. Otherwise, the hadronic phase is preferred over the quark-gluon 
phase. At the temperatures of interest, repulsive interactions reduce the contributions from 
the high mass part of the spectrum p3|, justifying our use of a finite number of resonance 



states in the EOS. It has been shown ||30| that within the different ways to include repulsion. 



a mean field approach, such as in Ref. gives the most realistic results. Thus, we will 



introduce a repulsive mean field potential V, which we couple only to the net baryon density. 
This is similar in spirit to the Walecka model . By including many of the resonance states, 
the main attractive interactions (akin to the scalar interactions of the Walecka model) are 
already taken into account. 

We assume that the repulsive potential energy density V is a function of pe of the form 

V(pb) = \^pI. (12) 
10 



where K is the mean field repulsion parameter. The partition function with a mean 
field interaction Eq. (|T^) is given by 



Z^^ = exp {-(3V[V{pb) - PbV'(pb)]} H ^xp [Zh{T, V, ^xf , /xs)J , (13) 

h 

where the effective baryon chemical potential 

< = /XB - V'(pb) = Pb - i^PB (14) 

describes the shift in the particle energy by Kp-^ due to the repulsive interaction. Using Eq. 
(^) with Eq. (|l^) leads to a self-consistency equation for the baryon density 



PB E / _ ^^5^ + _ ^g5,)/T] + 1 • (^^^ 

Once pb is solved, applying Eqs. ([10|) and (pJ]) to the mean field partition function, Eq. (|13|), 
gives the repulsion corrected energy density and pressure. The resulting EOS is labeled as 
EOS H and has one free parameter K. The ideal gas equation of state p(£:,pb) = £^/3, 
valid for massless non-interacting particles, is referred to as EOS I. We use this with three 
massless pions, to illustrate the effect of the number of hadronic degrees of freedom. 

So far, we have considered the equation of state using only the hadronic degrees of 
freedom. We will now construct an EOS having a phase transition to the QGP. As discussed 
earlier, the QGP in the critical region is highly nonperturbative and is best understood from 
QCD lattice simulations. For simplicity, the results of these simulations can be parametrized 
in terms of an ideal massless parton gas. In terms of temperature T and chemical potential 
PBi the thermodynamic densities are given by 

f 2 ^'^ 47r2 ^ ' 

e = 3p + 4:B (17) 

p^ = NfplT + -^pl, (18) 
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where /Xq = /Ub/3. B is the bag constant, and we use the number of flavors Nf = 2.5 in 
order to simulate efi^ects of the flnite strange quark mass. 



The phase boundary is determined by the pressure balance, phg = Pqgp, between the 
two phases at equilibrium. In the mixed phase, e and pe are calculated using the Maxwell 
construction. We define 

^QGP 

^(^'/'b) = ^HG + ^QGP 5 0<W<1 (19) 

as the volume fraction of the QGP in the mixed phase. 

The resulting EOS depends on two parameters, the mean field repulsion parameter K 
and the bag constant B. In this work, we give results for only two different choices of 
K and B, and call these parametrizations EOS A and EOS B. The two parameter sets 
are chosen to resemble the known features of the phase transition from lattice calculations 
26| , the difference encompassing the uncertainty of lattice results. These indicate a value 



Tc ~ 140 — 160 MeV for the transition temperature and Ae ^ 1 — 1.5 GeV/fm^ for the latent 
heat. Our EOS A, given by = 450 MeV fm^ and 5^/^ = 235 MeV, is similar to the one 
used in Ref. The resulting is 165 MeV, and Ae = 1.4 GeV/fm^ For EOS B, the 



parameters are K = 660 MeV fm^ and B^/^ = 200 MeV, resulting in = 140 MeV and 
Ae = 0.8 GeV/fm^. This parametrization should represent reasonable lower bounds for Tc, 
fic, and the latent heat. The parameters for the EOSs we use are summarized in Table |. 

The phase diagrams of equations of state A and B are shown in Fig. |l]. The difference 
in the latent heat can be seen if one plots the phase boundaries in the plane. This 
can be read off from Fig. 0, where the pressure is shown as a function of pb and e. For the 
EOS A, there is a large increase in the pressure, even in the mixed phase, with increasing 
baryon density, while for the EOS B the increase is much smaller. It would be interesting to 
obtain from hydrodynamical calculations constraints on the EOS at finite baryon density, 
since nothing is known so far from lattice calculations in this region. Presently, however, 
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the large uncertainties in the initial conditions preclude definitive conclusions to be drawn. 



IV. THE INITIAL CONDITIONS 

Hydrodynamical simulations at energies near or below 10 A GeV (BEVALAC and AGS) 
usually start with the colliding nuclei before the impact and include the initial compression 
and particle production. In the one-fiuid hydrodynamics, the nuclei fuse to a single fiuid, 
implying, at zero impact parameter, a complete stopping of equal size nuclei. At higher 
energies, as at the CERN-SPS, RHIC and LHC, this is not justified, and one must be able 
to incorporate nuclear transparency in the description. Instead of trying to describe the 
production and equilibration within hydrodynamics, we start the calculation from initial 
conditions which specify the hydrodynamic state of the system at time to- Initial conditions 
parametrize the production and equilibration dynamics. Note that to only plays the role 
of a bookkeeping device in the numerical calculations. Thermalization time does not enter 
explicitly in the parametrization, but, physically, it may be related to the initial longitudinal 
size, zq, at thermalization. In some cases, we find it useful to relate zq with to when defining 
the initial longitudinal velocity profile. 

Physically, the uncertainty in the initial conditions arises mainly from the lack of defini- 
tive knowledge about the nuclear stopping power and the time scale for equilibration. Two 
extreme scenarios are the full stopping model of Landau |]T| and the full scaling expansion 
model espoused by Bjorken |^|, in which Vz = z/t, beginning with equilibration. Although 



the precise energies at which these extreme cases are practically useful is unknown, the first 
case is expected to apply at moderately high energies and the second at ultrarelativistic 
energies. Since the present experiments fall between these limits, we have developed param- 
eterizations which span the range between these two extremes. We will try to incorporate 
some known features of hadron-hadron collisions in the parametrization, which are also 
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constrained to satisfy the conservation of energy and baryon number. 



We first consider the initial velocities. The collective four velocity is denoted by 
= 7(1, w), where v is the flow velocity vector. Since we consider zero impact param- 
eter cylindrically symmetric collisions only, we do not expect significant collective motion 
initially, and take the velocity in the radial direction at to to vanish, i.e., Vr{to, r) = 0. Note 
that for the Landau initial conditions, the initial longitudinal velocity also vanishes, i.e., 
Vz{to,T,z) = 0. Strictly speaking, the Bjorken model applies only in the infinite energy 
limit. In this case, the scaling ansatz for the four velocity is 



implying Vz = z/t. Invariance under longitudinal Lorentz boosts means that the thermody- 
namic quantities depend only on the longitudinal proper time r = — z'^, which equals 
the local time in the rest frame of any fluid element. 

At finite, albeit high collision energies, the longitudinal extent of the system is finite and 
the Bjorken scenario does not work properly in the fragmentation regions. Since we perform 
the numerical calculation using variables in the center-of-mass frame of the participating 
nucleons, we specify the initial condition at a fixed time to in this frame. We choose an 
initial velocity profile of the following form: 



In this form, should be regarded as a constant, which fixes the rapidity of produced matter 
at the edge zq, rather than the equilibration time. 

The reason for taking the rapidity y, instead of velocity Vz, to be proportional to z 
is purely practical. For numerical calculations, initial conditions must be smoothed and 
extended over the edge of the produced matter, initially at zq. The above parametrization 
leads to the Bjorken limit, Eq. pO|), in the inner part {z/to ^ 1), and extrapolates the 




(20) 





(21) 
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velocity smoothly to unity in the outer parts, where the thermodynamic densities approach 
zero. 



The energy density distribution in the Bjorken model is constant along contours of equal 



proper time ro(t, z) = \Jt'^ — . However, in the global frame, at fixed time t = to? it scales 
with z as [eq := e{to, 0)) |Q 



4/3 



7(^)'/' , (22) 



if an ideal gas EOS, e = 3p, is assumed. In the last term, 7(2;) = (1 — (^/to)^)^^^^ is the 
relativistic 7-factor for the Bjorken expansion velocity, Eq. 



We use Eq. ( p2D for our initial parametrization, despite the fact that it is strictly valid 
only for the velocity field Eq. ([20|), i.e., the Bjorken picture. One factor of 7 in this equation 
expresses the time dilatation effect for the moving cells and the factor 7^/^ the energy loss 
due to the work done by pressure in the expansion. Retaining this factor whenever the 
initial longitudinal velocity is non-vanishing is therefore physically well motivated. Then 
the remaining factor in e{tQ,z) represents the energy density in the local rest frame, and 
values at different z have the same physical interpretation. In the Bjorken picture, a rapidity 
plateau indicates constant local energy density. At fixed global initial time, this leads to 
increasing e{tQ, z) with increasing z. Since, at finite collision energy, the thermodynamic 
densities vanish when z > zq, this constraint will ensure that the energy density will have 



a maximum at finite z . This is in contrast to Ref. [18|, in which the authors employed a a 



constant energy density along z until the edges are reached. 

In production processes, distributions are always cut off smoothly when the phase space 
boundary is approached. Smoothing also helps to avoid oscillations in the numerical calcu- 
lations. We implement smoothing by multiplying the distributions with a Fermi function 

f{x,Xo,ax)= 777—. r-, f— r , (23) 

exp[(|x| - Xo)/a^J + 1 

15 



where a^, x = r,z is the diffuseness parameter in the radial and the longitudinal directions, 
respectively. The initial energy distribution is then given by 

e{z, r) = eo 7(z)^/^ f{z, zq, a,) f{r, tq, a,.) . (24) 

Since our code is symmetric in z in the center-of-mass frame of the participant nucleons, or 
equivalently, of the produced fireball, we cannot reproduce the asymmetric shape of rapidity 
distributions in S+Au collisions. In the central region, y = 0, we expect effects from the 
asymmetry on the expansion dynamics to be less important than those from the average 
longitudinal gradients, which are properly accounted for in the code. 

In the case of the longitudinal scaling expansion, the rapidity density of baryon number, 
dN^/dy does not change during the expansion. This allows us to relate the rapidity density 
to the initial baryon density through the equation 

T dy 



1 dNB dy , 1 dNE 

Pb[z) = —pyj — 1 — -r-: —r- T dz 



7(2;) dz 



(25) 



7TR'^^^^dz'-f{z) dy dz TTRpj.^^T dy 

where, for the scaling expansion, the expression in square brackets equals unity. In our case, 
expansion modifies the rapidity distributions and the problem is to find a reasonable choice 
for the initial dN^/dy. 

In nucleon-nucleon collisions, the leading particle effect in the nucleon distributions 
causes more nucleons to be present in the fragmentation regions than at central rapidity 
|3^,40|. In nuclear collisions, however, more stopping is present and it depends on the mass 



numbers of the colliding nuclei. Here, we assume that the initial rapidity density of baryons 
is flat, implying constant density at constant proper time in the Bjorken parametrization. 
At fixed global time, the baryon distribution as a function of z is then 

Psiz) = bo 7(2;) f{z, zo, /(r, tq, a^) , (26) 

where 60 is the central baryon density. 

Some of the initial parameters can be fixed by geometry. A natural choice for tq is the 
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sulphur radius, 3.65 fm. For the diffuseness parameter in a Woods-Saxon |^T[ parametriza- 
tion of the nuclear surface, we take the value = 0.5 fm. We set = 0.13 fm in order to 
have a^/ro ~ clz/zq. The number of participating nucleons is calculated from geometry, by 
assuming that at impact parameter zero the sulphur nucleus coalesces with a central tube 
of the gold nucleus having the diameter of sulphur. From the calculated value Stot = 104, 
we obtain the rapidity of the cm frame of the participating nucleons, ycm = 2.6, and the 
total cm energy, Etot = 970 GeV. This total energy should be considered to give the overall 
scale of the energy in the fireball, since small changes in the rapidity distributions of final 
particles in the fragmentation regions easily cause 10-20% changes in the total energy. We 
do not consider such deviations to be significant when fitting the pion spectra. Given Etot 
and -Btot; we can fix two of the remaining five parameters in the initial conditions. These 
are the central initial energy density Eq and baryon density 6o- 

The real ambiguity in the initial condition lies in the parameters zq and to- While in the 
Bjorken scenario, they are related through to = zq/c, i.e., the edge of matter is by definition 
on the light cone; this does not hold in a finite energy scenario, like the one introduced here. 
Now, the initial longitudinal velocity can deviate from the scaling behavior, f^^sc = z/t, and 
with our parametrization, Eq. (|21|), fixing to is equivalent to fixing i/q, the rapidity at Zq, as 



is done also in Refs. 115,18 . 



The values of Zq and to are chosen to fit the pion spectra. In Fig. ^, we show the 
initial profiles for the calculation with the EOS A. The main features are the maximum 
of the energy density at finite z instead of at 2; = 0, and a smooth behavior instead of a 
discontinuous one at the edge of matter. During the calculation, the velocity of the first 
vacuum cell is set equal to the velocity of the nearest cell containing matter. 

The profiles of the energy density for the other EOSs look similar. The parameters are 
listed in Table |I|. In calculating the total energy and the baryon number of the initial 
matter, only the volume with e{x) > Sf (£/ is the freeze-out energy density) is included in 
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the integration. The results are equal to the total energy and baryon number fluxes through 
the freeze-out surface. 



V. FREEZE-OUT 

We deflne the freeze-out surface as the space-time hypersurface cr^ of constant energy 
density. The value of £/ is chosen in such a way that the mean value of the temperature on 

is Tf ~ 140 MeV. This value for the temperature is chosen on the basis of calculations 
comparing the various mean free paths of the hadrons in the system with the size of the 



fireball 42-45 



The hadron spectra are calculated by assuming thermal momentum distributions and 
chemical equilibrium on the freeze-out surface. The invariant momentum distribution of a 



hadron h is then given by [46 



E ^ = ^ / da.jf , (27) 

(Pp {2-Kf I exp[(p^M'^-/i)/r/]±l ' ^ ' 

where the temperature Tf{x), the chemical potential fi{x), and the fluid flow four-velocity 
u'^{x) are determined on the surface from the hydrodynamic calculation. After calculating 
the momentum distribution for each hadron included in the EOS, Eq. the contributions 
from unstable resonance decays are added to the stable hadron spectra. We use the ap- 
proximations and decay kinematics described in Ref. |^. Finally, we integrate over the 
experimental acceptance in pt and y. 

VI. PHOTONS 

The thermal emission rate for photons can be shown to be directly proportional to the 
trace of the retarded photon self energy at flnite temperature . In the QGP, the imaginary 



parts of the lowest order contributions to the self energy correspond to tree level QCD 
Compton and annihilation diagrams, qq (77, q{q)g q{q)li respectively. Contributions 
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from tliese diagrams alone are infrared divergent. However, it has been shown by Braaten 
and Pisarski |^ that these long range effects are screened at finite temperature. The above 
mentioned diagrams were calculated, including the Braaten-Pisarski resummation, in Ref. 
||50| , |5T[] for zero baryon chemical potential. This result was extended in Ref. ||5^ to finite 
baryon density. We use the results of Ref. , since we have finite baryon density explicitly 



in our calculation. However, the influence of the term containing the chemical potential is 



small. For two quark flavors, the rate is 



where dR = dN/d^x. We use a temperature dependent running coupling constant 

" (33-2n;)ln(T/AT) ' ^^^^ 



where we take = 40 MeV |33| and, as in the EOS, nj = 2.5. 



For the hadron phase, we use the single photon production rates calculated in Ref. [^0 
These calculations were performed using a pseudoscalar-vector Lagrangian of the form 



C = - - Ip^^pA^- + Imjp^p^ - \f,,F^^ , (30) 

where and are complex pseudoscalar, photon and vector fields, respectively. Fur- 

thermore, = ~ ieA^ — iQpP^ is the covariant derivative, and Ppi, and F^^, are the 
vector meson and electromagnetic field strength tensors, respectively. Also, 5'p/47r = 2.9 as 
determined from the decay p nn. 

Parameterizations of the rates contributing to E^dR^'^ /d^p were given as a function of 



T and E^ in Ref. |]53| for the most important processes. These include the two scattering 
processes tt + tt ^ p + 7 and vr + p ^ tt + 7. The latter was only calculated for virtual pion 
exchange. It was shown in Ref. |B3 that 7 production in the vr-p channel is dominated by the 



ai decays. Thus, we include the process 7r + p-^ai^7r + 7by using the parametrization 
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given in Ref. ^EM. The decays p inT'y and u tt^'J, which occur during the hfetime of 



the fireball, were also included using the parametrization given in Ref. 



The total single photon spectrum is given by integrating over the total space-time volume, 
i.e., over all fluid cells with e > ej: 

E-— = / d x\w{e,pB) E— — {p-u,T,fiB) 

+ [l-w{e,ps)]E^{p.u,T)}, (31) 

where w is defined in Eq. (|19|). One should note that the rates are functions of p ■ m = p'^u^, 
the energy of the photons in the rest frame of the emitting fluid element. So far, the rates for 
the hadron phase do not contain processes involving baryons. Hence, there is no dependence 
on /iB. Possible contributions from baryons to the thermal photon yields were estimated to 
be small, even for Pb + Pb collisions at the CERN-SPS in Ref. |]55[ . 

VII. DILEPTONS 

We turn now to a discussion of dilepton production from both the quark-gluon plasma 
and the hadron gas. The dominant process in the QGP is the reaction qq l^l^, which 



was computed in lowest order for finite baryon chemical potential to be ||56 



dM^dPp 9 Svr 



2 



1 + 2ml /M^) V 1 - 4m2/M2 e"^/^ 



1^ {x- + exp[-(E + /i)/T]}[x+ + exp(-/x/T)] 

{x+ + exp[-(E + /i)/T]}[x_ + exp(-/i/r)] ' ^ ^ 

where x± = exp[-(^ ± |p|)/(2T)]. 

At low invariant masses, reactions at order 0{a'^as) become important [5?|. However, 



the mass region where these corrections are significant is dominated by the Dalitz decays of 
the final mesons [Q. We have not taken these higher order contributions into account in 
our hydrodynamical simulations. 
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For the hadron gas contribution to the dilepton rate, we use the results of Ref. ^9 . 
As in the calculation of the photon rates, the authors made use of the most general lowest 
order Lagrangian with vector and pseudoscalar mesons. The coupling of these fields to the 
electromagnetic current was computed under the assumption of Vector Meson Dominance. 
Prior to these calculations, the only rates computed were those due to tttt annihilation. 
The addition of other vector, pseudoscalar, and axial channels enhances the total rate by 
at least an order of magnitude away from the mass region of the light vector mesons. In 
doing so, however, one must not overcount channels already accounted for in the basic vrvr 



reactions [Q. Therefore, the vector meson decays were excluded |^T]] from the total rate 



given in Ref. ^9|]. We have also verified that the rates used in our calculations are in essential 



agreement with those calculated recently in Ref. using a spectral function approach. 



We adopt the following procedure in incorporating these rates in our hydrodynamic 
simulations. Our starting point is the parametrization of the total dilepton production rate 
dR^°^ /dM'^{M,T) as a function of temperature T and invariant mass = = {pi +^2)^, 
as in Ref. |5^. Here, p, pi, and p2 are the four-momenta of the pair and the single leptons, 
respectively. 

In order to apply the kinematical cuts, we need the momentum distribution of the pair. 
To obtain this, we use the relation 



dR ^ 1 E IT dR . . { oo\ 

dMHypTdpT ^ 2MTKi(M/T) ~dM^^ ' ^ ^ 

where E = p^u^^ is the pair energy and Ki is a modified Bessel function. This relation is 
valid for reactions in which the final state contains only the lepton pair, which gives the 
dominant contribution in the higher mass region. However, the relation is not valid for the 
decays, h —>■ h'e^e~, which result in small mass pairs. The region in which Eq. ( ^31) is a 
good approximation was given to be above 300 MeV in invariant mass in MM. We have 



verified that the approximation Eq. (pSQ adequately reproduces the rates in Ref. |^ above 



an invariant mass of 400 MeV. Below 400 MeV, differences may exceed a factor of 2 and 
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above in some regions of the phase space. However, in this region the spectrum is dominated 
by Dahtz decays of mesons after freeze-out, which are accurately treated separately. 



The rates, Eq. (^) and Eq. (|33D , are integrated over the space-time volume of the 
fireball, as in Eq. ( ^1]) for the photons. Here also, we include baryon contributions from only 
the QGP phase. It remains to be seen whether baryons can contribute significantly to the 
dilepton yields (see for example, Ref. for initial estimates). 

The measured dilepton spectra also contain pairs from decays after freeze-out. We con- 
sider this as a background contribution. When comparing our results with the measurement 



of the CERES-coUaboration [14|, we calculate this background from our hydrodynamical 



simulation, instead of using the background estimated by the CERES-collaboration [p!4[ . 
We combine the freeze-out momentum distributions of 7r°, rj, p°, u, rj', and as given by 
Eq. ( p7D with the distributions from the decays of higher lying resonances. 

First, we consider the vector meson decays into electron pairs. The hydrodynamic cal- 
culation is done assuming a fixed mass and zero width for the resonance states. In order 
to calculate the dilepton spectra, we have to take into account at least the width of the p° 
meson. We assume that the mass distribution of the resonance is of the Breit-Wigner form: 

dM^d^p 71 {M^ -m^y + rl,{M)ml. d^p ' ^ ' 

where mv is the mean mass of the vector mesons. The normalization ^ is determined in such 
a way that the total yield has the value calculated at the freeze-out. We use the following 



mass dependence of the total p width ||47[| : 



with M and in units of GeV. For the other vector mesons [uj, </>), constant widths taken 



from experiment |66] are used. This is justified, since the experimental resolution of CERES 



is much wider than these widths. In addition to the cuts 2m7r {2m k) for p° (0) at threshold, 
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we apply a cut Sm^r for the u mass distribution. The electro-magnetic decay width in the 
Vector Meson Dominance model is proportional to the pair mass [Q. Hence, we use the 
parametrization 

M 



my 



(36) 



with the values of my and Ty^^+g. taken from Ref . |^ . Thus, the electron pair distribution 
is given by 



(37) 



Next, we consider the Dalitz-decays of the mesons. We use the differential decay width 



dr 



dM^ 



3^(l+2</M^) 

2 \ 2 



1 + 



1 - Aml/M^ 



(mt 



m 



2\2 



(38) 



If the particle j is a photon 7, the normalization to Fj^^^ results in an additional factor of 

2 on the r.h.s. of Eq. (pS]). For tt^ 7e^e~, we take a linear approximation for the form 

factor, F(M2) = 1 + 5.5 GeV^^ x , from Ref. For 77 ^ 76+6" and cu 7r°e+e-, the 

dipole approximation, ^(M^) = (1 - MVA^)~\ with = 720 MeV and = 650 MeV, 

respectively, is used. For rj' — > ■je~^e~ , the Vector Meson Dominance form factor 

,4 



|F(M 



2m2 



is used. The pair distribution is then calculated from 

»^je (M) (M) 



(39) 



(40) 



dw^d^p r?^p rfM2 ' ' d^p 

where Ftol* is taken from Ref. ||6^ and EdN^^^'^^ /d^p is an isotropic electron pair distribution 
resulting from the decay of an unpolarized meson i into a pair of mass M and a particle j', 
as given in Ref. |^^. The incorporation of the experimental cuts and resolution is described 
in the appendix. 
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VIII. RESULTS 



In this section, we discuss the results of our attempt to describe consistently the hadron 
and electromagnetic spectra in S+Au collisions using the hydrodynamic approach. We 
fix the parameters for the initial conditions from the hadron spectra. In general, good 
fits are obtained only for the pion distributions — the fits to heavier hadron spectra are 
significantly worse. These are discussed in detail below. We apply the temperature and 
velocity distributions corresponding to these fits to predict the photon and dilepton spectra 
for each EOS. Our results may be summarized as follows. The existing photon data are 
not sufficiently precise to exclude any EOS except the EOS for an ideal, massless pion 
gas. Hopefully, future experiments will allow for better discrimination. With regard to 
the CERES dilepton data, none of the EOSs considered, in conjunction with the standard 
leading order dilepton rates, succeed in reproducing the observed excess of dileptons below 
the p peak. 



A. Hadron Spectra 

We begin our discussion with the results for the hadron spectra in S+Au collisions. In 
Fig. ^ we show the rapidity spectra for several hadrons obtained from calculations using 
three different EOSs. In choosing the parameters in the calculation, our main emphasis has 
been to reproduce the negative particle spectrum. In the case of the ideal pion gas EOS I, 
we assume that the negatives consist of negative massless pions only. For the EOSs A and 
B, both of which exhibit a phase transition, the initial conditions are quite similar, but for 
the two purely hadronic EOSs H and I, they differ considerably (c.f. Table IH). 

Let us first note the features that are common to the three different EOSs. For the 
negatively charged particles, which are mostly pions, the calculations are in good agreement 
with the data from NA35 0. This tells us that, by a suitable choice of the initial conditions, 
the flow of energy density across the freeze-out surface can be reproduced reasonably well 
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with all the EOSs considered. The three pion states carry roughly half of the total energy. 
The other half is carried mainly by the baryons. 



The experimental net proton p — p distribution in Fig. ^ shows a large asymmetry with 
respect to the center of momentum. (The two data points with y < 3 are identified proton 
data from a different experiment NA44 with S + Pb. It is instructive to view these 
points together with the net proton distribution from NA35 0, because in this region the 
anti-proton yield is very small. While NA35 has practically a full pr-acceptance up to 2 
GeV, the NA44 data were extrapolated to the full pt region.) There is a large difference in 
the target and projectile size, and we are presently unable to take this into account in our 
code with longitudinal mirror symmetry. However, we do not want to neglect the baryons 
totally, since they carry a significant part of the total energy; their influence is felt through 
the EOS during the space-time evolution. We reproduce the right total number of baryons in 
our calculation. Also, the baryon rapidity density in the central rapidity region {ycm = 2.6) 
is roughly reproduced. This should ensure a reasonable description of the expansion in 
the central region, where our calculations of the electromagnetic yields are compared with 
measurements. 

The Kg distribution comes out too large in comparison with experiment. This may be 
attributed to the fact that we have assumed full chemical equilibrium. A detailed chemical 
analysis of strange particle spectra in the S + Au collisions in Ref. ||3^ shows a strangeness 



suppression of 0.7 relative to full chemical equilibrium. Multiplying the i^°-spectra and the 
other strange particle spectra by this factor would improve the fit. The abundance of 
mesons is also slightly influenced by the baryon distributions through the requirement of 
strangeness neutrality. The and have to compensate the net negative strangeness 
of the A baryons. In the region 3 < y < A, the overestimation of the A's leads to a small 
additional contribution to the K^^s. 

The A and A rapidity distributions are influenced both by the strangeness and the baryon 

25 



number chemical potentials. From the rapidity distribution of the protons, we see that the 
baryon number, or equivalently the baryon chemical potential, is too large in the forward 
rapidity region, because of the longitudinal symmetry in our calculation. The excess in the 
A spectrum in the region 3 < ?/ < 4 is similar to that for protons. In order to answer 
the question of how well strangeness chemical equilibrium holds for the strange baryons, 
we should first reduce /ib to reproduce the proton distribution in this rapidity range. It 
seems that this would lead to a reasonable agreement for the A's, assuming that fis would 
be unchanged. However, a reduction in /i^ would not be sufficient to increase the A yield 
to the experimental value at rapidities around 3. 

We may summarize our analysis of the strange hadrons by stating that, for the freeze-out 
parameters that fit the negative particle yields, there is agreement in the yield of A's relative 
to that of protons, but the A's are underestimated and the Kq^s are slightly overestimated. 
An analysis of S+S collisions at the same freeze-out temperature, Tf = 140 MeV, has led to 



the same conclusion 18,68 



The particle yields are determined in our model by the assumption of chemical equilib- 
rium until freeze-out. This is poorly justified in the late dilute phase, since the strangeness 
changing cross sections are small. Therefore, a study of the chemical behavior in the late 
phase deserves further investigation on its own right, but is outside the scope of this inves- 
tigation. However, looking at the pion distributions, we think that the bulk behavior of the 
longitudinal expansion is well reproduced by our model. 

The transverse momentum distributions are shown in Fig. |^. The overall agreement is 
quite good with the exception of the EOS I. (We will return to discuss this case later.) The 
m-r-spectra of vr'^s are well fit up to 2 GeV. For larger transverse momenta, the three EOSs 
can also be distinguished. If one uses the same initial conditions, then a stiffer EOS would 
produce more transverse fiow. However, to fit the longitudinal spectra, we enhanced the 
initial energy density for the EOS B and reduced it for the EOS H in comparison to that 
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for the EOS A (c.f. Table H). The initial temperatures are roughly the same in each case, 
and lead to the similar final py-slopes. 



Although the best agreement is achieved using the EOS H, fits for the other EOSs can 
be improved by fine-tuning the initial conditions. Further, it is doubtful that particles with 
very large pt follow a hydrodynamical behavior. Instead, they might originate from semi- 
hard processes. We therefore conclude that, with the exception of the EOS I, we cannot 
rule out any of the remaining three EOSs on the basis of fits to the pr-spectra. 

The ideal pion gas EOS I produces too fiat a slope. The initial conditions in this case 
correspond to the same total energy as in the other three cases (see Table H). This was 
necessary in order to get the rapidity spectra correct in the central rapidity region. Since 
there are no baryons to share the energy, the surplus energy is converted to transverse kinetic 
energy. In short, for the EOS I, we could not find initial conditions which reproduce both 
the rapidity and p^-spectra simultaneously (see also Ref. [j70|). Therefore, the ideal pion gas 
EOS I can be ruled out, on the basis of hadronic data, as being too stiff. 

The absolute yields of the other particles in Fig. ^ depend on the assumption of chemical 
equilibrium, the inadequacy of which we have discussed above. This is also corroborated 
by the the anti-proton yield, which is also underestimated, like that of the A's. We see, 
however, that the slopes for the different particles are reasonably reproduced by all three 
EOSs, supporting the picture of collective transverse flow present at SPS energies ll^ . 



At low Pt, especially for the A's, discrepancies between calculations and data persist. 
The calculation is expected to overshoot the data on the basis of the excess seen in the 
rapidity spectrum in the 3 < y < A region. The larger relative weight from this region in the 
calculation as compared to the data might also be a reason for the difference in the slopes. 



The experimental r^/vr -ratio for central events |Tl| has large error bars. It is reproduced 
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reasonably well at large pt, but our calculation underestimates it at low px- In our model, 
the 4:71 integrated r^/vr'^-ratio depends mainly on the freeze-out temperature. In Table |T|, we 
compare the total multiplicities with p + p data on meson production |p9|. Our r^/yr^-ratio 



is similar to that in proton collisions and the minimum bias data in S + Au collisions ||T2|. 



To summarize the discussion of the hadron spectra, we see that with the exception of the 
ideal pion gas EOS, the initial conditions can be adjusted to reproduce the gross behavior 
of the hadron (negatives) spectra for a wide class of EOSs. Details of the spectra depend on 
the assumption of chemical equilibrium, which turns out to be poorly justified at freeze-out 
temperatures of Ty ~ 140 MeV. Chemical equilibrium between particles and resonances 
which have a large cross section may still hold. 



B. Photon Spectra 

The fundamental difference between the hadron spectra and the photon spectra is the 
fact that photons are emitted from collisions of charged particles during the entire expansion 
stage. We have calculated the photon p^-spectra using the same simulation from which the 
hadron spectra were obtained. We integrate the photons over the four-dimensional space- 
time volume, Eq. (|3TD, which is bounded by the three-dimensional freeze-out surface used for 
the calculation of the hadron spectra. The calculations are compared with the upper bound 
of the WA80 collaboration on the direct photon spectrum [|T^]. Experimentally, the direct 
photon yield is obtained by subtracting photons from the decays of mesons and baryons. 

Since the photon yield depends on the properties of the system as it expands, we first 
discuss how the space-time behavior is affected by the different EOSs. In Fig. |^, contours 
of constant energy density in the ^t-plane at r = are shown, for all four EOSs. The 
boundaries of the mixed phase are indicated as solid lines. The freeze-out times at the 
center z = are tf = 9.4, 7.9, 7.9, and 7.5 fm/c for EOSs A, B, H, and I, respectively. 
At r = 0, transverse expansion is absent essentially up to tf. Similar freeze-out times for 
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different EOSs indicate that the longitudinal expansion in the central region is dominated 
by the initial velocity gradient. At large 2;, the QGP equations of state produce a long- living 
tail. This is the result of a weaker longitudinal acceleration, due to the smaller pressure 
gradient in the mixed phase as compared with that in a calculation employing the hadron 
gas EOS without phase transition. 

The space-time volume of the mixed phase in the EOS B is somewhat larger than that for 
the EOS A, despite the fact that it has a smaller latent heat, which allows, per unit volume, 
a faster conversion of the mixed phase into the hadron gas. However, for the EOS B, the 
mixed phase is reached at a later time and with so much larger volume that the conversion 
takes more time and produces a bigger final volume than that for the EOS A. 

In Fig. we show the direct photon spectra for the EOS A. From the individual 
contributions from the different phases, it is clear that hadronic processes dominate the 
production of single photons in this case. In Ref. [0], it was shown that at the same 
temperature, emission rates per unit volume are roughly the same in a QGP and a hadron 



gas. Inclusion of the ai mesons enhances photon production for energies > 1 GeV 
leading, at T^, to a considerably larger photon production rate in the hadron gas than in the 
QGP. Also, the high temperature phase above lasts for a short time in this case, and the 
produced matter spends most of its lifetime in the hadron gas phase. These two features 
lead to the dominance of photon emission from the hadronic phase. The higher temperature 
of the initial QGP phase shows up as a flatter px slope in the plasma contribution, but even 
at the largest pt values, it is clearly below the hadron gas contribution. The total yield is 
in agreement with the upper bound provided by the data. 

In Fig. the individual contributions are shown for the EOS B. In this case, the 
hadrons from the mixed phase dominate photon production. The QGP yields are similar 
to those with the EOS A. At large pt, the QGP contribution is as large as the hadron 
contribution from the mixed phase. There is contribution from the pure hadron gas phase, 
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because the transition from the mixed phase to the hadron phase and freeze-out takes place 
simultaneously at the same temperature. 



Figure || shows the total single photon spectra for the four different EOSs. The situation 
is qualitatively similar to that for the pr-spectrum of vr'^'s in Fig. ^. However, quantitative 
differences exist. A close inspection reveals that not only the results for the EOS I, but also 
for the other three EOSs, differ from one another. The total yield with the EOS B lies a 
factor 2-3 below that of the EOS A, whereas the EOS H produces 2-10 times more photons 
depending on the p-r-region. The yield for the EOS I lies orders of magnitude above. 

In principle, the dependence of the photon production on Tc could be used to determine 
Tc from the data. However, one can see from Fig. |^ that the present upper bound on the 
single photon yield rules out only the ideal pion gas EOS 1. (Recall that this EOS is also 
ruled out by the transverse momentum spectrum of negative hadrons.) From our results 
of the photon spectra, we cannot claim evidence for a phase transition, in contrast to the 
claims in Refs. [p^-|72[. The main point here is that, if a reasonable amount of degrees of 
freedom are taken into account in a hadron gas, the increase of temperature with energy 
density is reduced. This is clearly seen by comparing results for the EOS H with those for 
the EOS I. The large photon yield in the case of the EOS I is due to the very large initial 
temperature, Tj = 400 MeV, as compared to Tj = 250 MeV for the EOS H. 



In the extreme case of the Hagedorn bootstrap model |31|] , we have a limiting temperature 
~ 150 MeV in a hadron gas, with arbitrarily high energy densities. For limited energy 
densities, such an EOS would lead to a small photon production. Therefore, the results do 
not readily attest to the existence or the absence of a phase transition. With the present 
precision, one can only rule out high initial temperatures. However, since the rates from the 
hadron gas and the QGP differ, the correlation between the total yield and the slope of the 
transverse momentum distribution will differ for the purely hadronic EOS and an EOS with 
a phase transition. For this reason, improving the experimental upper limit and measuring 
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the yield of direct photons is very important. 



In discussing the dependence of the photon production on the transition temperature 
Tc, one should keep in mind not only the limitations of hydrodynamics, but also the uncer- 
tainties in the rate calculations both in the QGP and a hadron gas near Tc. Our discussion 
above is based on the considerable difference in the rates at the same temperature between 
perturbative QCD results |Q and leading order estimates in a hadron gas . 



1. Comparison with other works 

Finally, we want to compare our results with other calculations of photon yields in the S 
Au collisions, which have used hydrodynamics to describe the evolution of the produced 



matter. These earlier calculations were compared with the preliminary data of WA80 |]10 
The new analysis by WA80 gives upper bounds, which are compatible with the preliminary 
data. 



First, we confirm the result of Ref. that the absolute normalization of the photon 
yield is sensitive to Tc and to the EOS. The calculation in Ref. was done using a one- 
dimensional Bjorken-like hydrodynamical scenario, whereas we employ a three-dimensional 
expansion. The small differences in the photon yield from our results can be attributed to the 



differences in the hydrodynamical solutions. In Ref. [0, consistency with the preliminary 
data of WA80 was achieved with Tc ~ 200 MeV, while lower critical temperatures of around 
150 MeV led to an underprediction of the preliminary WA80 data. Therefore, the possibility 
of a long-lived mixed phase, of duration 30-40 fm/c, was suggested there. For the EOSs used 
in our three-dimensional calculations, Tc lies in the range 140-165 MeV and the duration of 
the mixed phase does not exceed 10 fm/c. Since the present data provides an upper bound 
only, a long-lived scenario is not necessary with the lower values of Tc which we have been 
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using. 



Our results roughly agree with the investigation in Ref. |^ as well. However, Arbex et 



al. |23| only investigated an EOS with Tc of 200 MeV. Since, as noted above, the absolute 
normalization of the photon spectra is sensitive to the critical temperature, their results 
agreed with the preliminary data of WA80. 



In Refs. |]70| and |7l[], an agreement with the preliminary data of WA80 was achieved using 



an EOS with Tc around 160 MeV, while our calculation with EOS A would underpredict this 



preliminary data by roughly a factor of 4. In Ref. , a rather low freeze-out temperature, 
Tf = 100 MeV, was chosen. We have checked that the photon yield increases by a factor 
of ~ 1.5 when the calculation is extended from a freeze-out temperature Tf = 140 MeV 
to Tf = 100 MeV. In Ref. |]7l|, the neglect of baryons and the smaller number of mesons 
in the EOS leads to considerably longer lifetimes for the mixed and hadron phases. In 
addition, our three-dimensional calculations lead to more rapid cooling than that obtained 
using the scaling expansion, even when transverse expansion is included. A combination of 



these effects can explain the larger yield in Ref. . 



C. Dielectrons 



The other electromagnetic signal, measured in the S-I-Au collisions by the CERES col- 



laboration ll^], is the dielectron mass distribution. In these measurements, the dilepton 
background from the decays of final mesons is not subtracted, because the amounts of dif- 
ferent mesons are not precisely known. Thus, we have two parts in the final dilepton yield. 
First, the emission during the lifetime of the fireball and, second, the electromagnetic decays 
of hadrons after the decoupling. The latter contribution is shown in Fig. |^ for the EOS A. 
This background is similar for the other EOSs, since the final hadron spectra are reproduced 
in each case by tuning the initial conditions, as discussed above. The calculation of the decay 
dilepton spectrum is described in Sec. VII and the appendix. All hadrons which produce 
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lepton pairs are considered. Both the thermal contribution and the contributions from the 
decays from heavier hadrons are included. 



Our calculated background in Fig. |^ is basically in agreement with the estimated back- 
ground in Ref. [14]. However, there are some differences in the total yields of the decaying 
mesons, since in our calculation the meson multiplicities are given by the calculated freeze- 
out conditions, mainly the assumed temperature, the effect of the baryon chemical potential 
being small. On the other hand, the CERES collaboration used the meson-to-7r° ratios from 
p+p collisions [B^ to fix the multiplicities of mesons from the measured 7r° spectrum [llO 



In Table [ITT| , we present the meson yields normalized to the vr" yield. Our rj multiplicity is 
somewhat higher than in p + p collisions, while that of the p° is slightly smaller. For the lo 
meson, the ratio is 3/4. The main differences in the input yields are those for the 77' and 
(p. The contributions from the 77' are negligible compared to those from the other dilepton 
sources. Thus, the only important difference with the CERES background is a factor of 4 
in the 0-mass region; however, this is still in agreement with the data. 



In Fig. p!0|a, we show the dileptons radiated during the lifetime of the fireball. These 
results are folded with the CERES cuts and the CERES resolution. Here, we see the 
same systematics as in the case of photons; the hadronic contribution dominates the yield. 
The largest contribution comes from the tt-tt annihilation to dielectrons via the p-form 
factor (Vector Dominance); this is seen as a peak at the p mass. As for the photons, the 
contribution from the pure hadron gas and the hadronic contribution in the mixed phase 
are equally important. 



The sum of the background and the thermal emission is shown in Fig. [TpD. We first note 
that thermal emission roughly fills the gap between the background and the data around 
and above the p mass, even though the calculated results tend to lie at the lower bound of 
the errors, especially for the EOS B. Note that the systematic and the statistical errors are 
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shown separately 







In the mass region between 200 and 600 MeV, our calculated results lie clearly below 
the data. While the calculation has a dip in this region between the contributions from the 
Dalitz pairs and the vector mesons, the data is flat and smooth. There have been several 
suggestions for the origin of the excess over the expected sources in this invariant mass region 
^,75-78|. An interesting possibility for the explanation is a shift of the vector meson masses 



associated with the expected restoration of chiral symmetry as the transition temperature Tc 
is approached in the hadron gas. We have not tried to include this in our calculation. With 
an appropriate parametrization of the temperature and/or density dependence of the hadron 
masses, we most likely would be able to reproduce the data, but a consistent treatment in a 
hydrodynamical approach would require the incorporation of density dependent masses to 
also calculate the EOS and the decay rates. This is beyond the scope of the present work, 
and its proper implementation requires a major effort, which will be taken up separately. 

Our conclusions for the EOS from the dilepton calculations are similar to those for 
photons. The absolute yields are sensitive to and the results for the EOS B with Tc = 140 
MeV are below the data at all mass values. However, the differences are not as pronounced 
as for the photons, since the contributions from the decays of final hadrons have not been 
subtracted. 



1. Comparison with other works 

The excess production of low-mass dilepton pairs in S + Au collisions was recently 
addressed in Ref. using a one- dimensional Bjorken-like expansion. In the invariant mass 
range 0.2 GeV < M < 0.6 GeV, where the excess over expected sources is most evident, 
the data were underpredicted by about a factor of 4. Due to the high initial temperature 
(Tj = 380 MeV) required to obtain similar results in a purely hadronic scenario, the case in 
which a QGP is admitted (T^ = 198.7 MeV, = 160 MeV, Tf = 140 MeV, and a mixed 
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phase duration ~ 13 fm/c) was favored in this work. 

Where the deviations from the data are largest, our three-dimensional calculations (with 
the same standard rates as above) underpredict the data by a much larger factor (of about 



8-10). Several sources for the differences from Ref. may be cited. These include the use 
of a more realistic EOS, a shorter duration of the mixed phase, and features specific to the 
three-dimensional flow of the matter emitting the lepton pairs. Clearly, a combination of 
these effects results in dilepton yields that are lower than in the case of a one-dimensional 
Bjorken evolution, even in the case where an EOS admitting a phase transition to the QGP 
is used. We also calculate the background contributions from the calculated hadron spectra. 



whereas the CERES background is used in Ref. [34 



A comparison of our hydrodynamical results with the alternative sequential scattering 
models (also termed as cascade or transport models) depends on the extent to which thermal- 
ization is achieved in the latter approach. Specific medium modifications of the vector meson 
properties, in particular a decrease in their masses, have been found to yield a satisfactory 
description of the CERES data | ]75|p5[ . Whether a similar approach can be satisfactorily 
adopted in hydrodynamical simulations is a challenging future task. 



IX. CONCLUSION 

The aim of this work was to establish the extent to which one can constrain the EOS 
from the experimental data for S + Au collisions from a simultaneous description of the 
hadron and electromagnetic spectra using hydrodynamics. We have shown that, in general, 
the influence of the EOS on hadronic spectra can be counterbalanced by choosing different 
initial conditions. A simultaneous calculation of the electromagnetic signals can, in principle, 
distinguish between the different EOSs. However, the present experimental resolution allows 
us to rule out only extreme cases, such as the ideal pion gas EOS with only a few degrees 
of freedom. Also, the dilepton yield for the QGP equation of state with Tc ^ 140 MeV, the 
EOS B, tends to fall below the data in the vector meson mass region, indicating an effective 
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lower bound of 140 MeV for Tc, if the transition exists. 

The constraint that can be drawn from the single photon data is that the initial temper- 
ature cannot be too high. The present data rules out temperatures above 250 MeV. This 
limit on the initial temperature can be achieved only if a large number of degrees of freedom 
is involved, be it in the form of quarks and gluons, or in the form of a large enough number 
of hadrons. However, if the data can be improved, the two cases can be distinguished, since 
the total emission from a hadron gas is larger than that from the QGP. In the total yield, 
the difference between a pure hadronic EOS and an EOS with a phase transition increases 
with decreasing Tc. 

The behavior of the dilepton spectrum in the mass region between 200 and 600 MeV 
shows that the description of the hot and dense strongly interacting matter near Tc in terms 
of the free-space parameters is not adequate. With our present hydrodynamical approach, 
the dilepton spectrum can be explained only in the mass region of the vector mesons. The 
large experimental dilepton yield below the p mass may indicate medium modifications of 
the particle properties. These effects can be included in hydrodynamical calculations, but 
for a consistent calculation, the EOS must be modified accordingly. 
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APPENDIX 



CUTS AND RESOLUTION FOR COMPARISON WITH CERES 

The kinematic cuts of the CERES experiment are cuts in the momenta of the electron and 
the positron We incorporate them in the following way. We take dN / {dM'^ dyprdpr) 



of the pair calculated either from the background contribution (see Eq. ^ and Eq. ^0]) or 
that resulting from the fireball during its lifetime (see Eq. (^2]) and Eq. (|33D). Going to the 
pair rest frame, we assume an isotropic momentum distribution. Thus, the single electron 
momentum distribution dN/ (dM'^dyipiTdpiT) , when the pair mass is M, is given by 

y+ (mj)^ 

M f ^ 1 ^ 2 at cut^ 

-TTTT-i 1 = . dy / dniT 0{p2T - Pt ) 

dM^dyiPirdpiT 2nJ{M^ - Ami) ~ 2 

X ^^(Cax - 1/2 - Vcr.) 9{y2 + - Cn) ^(^t " Q"""') 

^ 1 dN ^^^^ 

\IptP\t - [^72 - mrmiT cosh(?/i - y)f dM^dyprdpr ' 

where y, yi are the rapidities in the fireball rest frame. Further, 



y^ = yi± sinh-i[J(M2 - 4ml)/{2miT)] 



m. 



^ M MrriiT cosh(?/ — ± pixyM'^ — Ami ~ 4m^2. sinh^(|/ — yi) 



(42) 



2 m\rp smh^ {y — yi) + ml 

The opening angle of the electrons in the laboratory system, d^l2i neglecting the electron 
mass, is given by 

cos«t) = 1 - , „^ (43) 

^ ' 2 pi l^b p lab V i 



with \pi['''° = PiT\Jl + sinh^(?/i + j/cm). 

The final spectrum of the pairs in the approximation r^gi = ?/e± is then 

dA^™* 2M f f dN . . 

dMd-rj = ^J ^^'^ J dM^dy,p,rdp,T ' ^^^^ 

„cut ^cut 
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In order to compare with the CERES experiment, we finally have to fold the calculated 
results with the detector resolution. We use a Gaussian folding 

-(M) = / t/M'-^— — exp -i-— — f -7TT:riM ) ' (45) 



dMdrj ' ' J y/2na{M') \ 2a'^{M') J dMdrj 

with the width o"(M) taken to be the detector resolution kindly provided to us by the CERES 
collaboration llTii. 
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TABLE |. Some physical features characterizing various model EOSs. 





EOS A 


EOS B 


EOS H 


EOS I 


Hadronic 

Degrees of freedom 


p, n, r]', 0, A, S 
A,ai,S,E(1385) 


p, n, r]', 0, A, S 
A,ai,S,E(1385) 


TC,K,r],p,uj,K* 
p,n,r/',0,A,S 
A,ai,S,S(1385) 


vr 


Number of 

QGP degrees of freedom 


31 


31 






Mean field repulsion 
K [fm^ MeV] 


450 


660 


450 





Bag constant 
51/^ [MeV] 


235 


200 






Te [MeV] 


165 


140 


00 


00 


/iBc [MeV] 


1770 


1290 


00 


00 


Ae [GeV/fm3] 


1.4 


0.8 
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TABLE m Variables used for the parametrization of the initial and freeze-out conditions 
for the different EOSs. The first part contains the free parameters, the second part some 
deduced quantities, and the third part the freeze-out energy density, with the corresponding 
temperature and baryon chemical potential averaged over the freeze-out hypersurface. The 
symbols are explained in the text. 





EOS A 


EOS B 


EOS H 


EOS I 


ro [fm] 


3.65 


3.65 


3.65 


3.65 


ttr [fm] 


0.50 


0.50 


0.50 


0.50 


Zo [fm] 


0.80 


0.80 


1.2 


1.97 


az [fm] 


0.13 


0.13 


0.13 


0.13 


eo[GeV/fm3] 


7.0 


8.0 


5.0 


3.3 


bo [fm-3] 


1.10 


1.16 


0.75 





to' [c/fm] 


1.06 


0.95 


0.71 


0.41 


viz = Zo) 


0.85 


0.76 


0.85 


0.80 


Etot [GeV] 


914 


923 


904 


994 


-Btot 


105 


104 


104 





'S'tot 


3880 


4230 


3650 


3035 


S/B 


37 


41 


35 




T,{z = 0) [MeV] 


238 


249 


248 


400 


tf [fm/c] 


9.4 


7.9 


7.9 


7.5 


£dcc[GeV/fm3] 


0.165 


0.15 


0.165 


0.05 


(Tf) [MeV] 


142 


141 


142 


140 


(/iB/) [MeV] 


230 


210 


235 
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TABLE |T|. Total multiplicities used to estimate the background to the dielectron spec- 
trum. The total yields are normalized to the 7r° yield. The EOS A is used in the calculation. 
For comparison, the total production cross section in p + p collisions at 400 GeV incident 



momentum are shown in the first column (the data are taken from NA27 p9| ) 





Hydro 


NA27 




total 


relative 


total cross 


relative 




multiplicity 


to TfO 


section [mb] 


to 7r° 




212 


1.000 


127.2 


1.000 


V 


19.4 


0.092 


9.8 


0.077 


P' 


18.6 


0.088 


12.6 


0.099 


UJ 


15.8 


0.074 


12.8 


0.101 




2.0 


0.009 









4.1 


0.019 


0.62 


0.005 
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FIGURE CAPTIONS 



Fig. |T|. Phase diagrams for two equations of state with a first order phase transition. The 
meaning of the parameters is explained in the text. For the EOS A, K = 450 MeV fm^, 
B^/^ = 235 MeV, and for the EOS B, K = 660 MeV fm^, B^/^ = 200 MeV. 

Fig. ||. Pressure P as a function of baryon pe and energy densities e for the EOS A (a) 
and the EOS B (b). 

Fig. |. 

Initial conditions for S + Au collisions at 200 A GeV/c for the case of the EOS A. The 
different panels show (a) the rapidity y, (b) the longitudinal velocity Vz, (c) the energy 
density e, and (d) the baryon density pe as a function of the longitudinal coordinate z; (e) 
the baryon density pe as a function of rapidity y; (f) the energy density £ as a function of 
radius r. 

Fig. ^ Rapidity distributions for several hadrons compared with a symmetric hydrody- 
namical calculation. The equations of state are EOS A (solid line), EOS B (dotted line), 
hadronic EOS H (short dashed line), and ideal pion gas EOS I (long dashed line). The data 
are measured by the following groups: Negative charged particle data are from NA35 0, 
the strange hadron (A, A, K^) data from NA35 the net proton {p — p) data, with y > 3 
from NA35 and the two points of proton data with y < 3 are from NA44 

Fig. ^. Transverse momentum distributions for several hadrons compared with a symmetric 
hydrodynamical calculation. Curve designations are as in Fig. ^. The data are measured 
by the following groups: vr'^ data for central S + Au collisions from WA80 [|l^], the strange 
hadron (A, A, K^) data [Q, and the anti-protons are from NA35 |^, and the tj/tc^ ratio for 
central S + Au collisions from WA80 ||11|| . 

Fig. ||. Contours of constant energy density in the zt-p\a.iae for calculations with the EOS 
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A (a), EOS B (b), EOS H (c), and EOS I (d). Contours counted outward correspond to 
energy densities of 4.0, 2.0, 0.9, 0.4, 0.15 GeV/fm^ for (a) through (c), and to 2.0, 0.9, 0.4, 
0.15, 0.05 GeV/fm^ for (d). The sohd hues indicate the transition from the QGP to the 
mixed phase and from the mixed phase to the hadronic phase. The dashed hne corresponds 
to freeze-out. In (b), freeze-out occurs at the transition from the mixed to the hadronic 
phase. Note the different scale in (b). 

Fig. 0. Single photon py-spectra compared with the upper bound of WA80 data |]13[; 
(a) EOS A and (b) EOS B. The different contributions shown are: HG hadron gas, M(HG) 
hadronic part of the mixed phase, M(QGP) QGP part of the mixed phase, QGP, and the 
total spectrum (solid line). 

Fig. ^. Total spectrum of single photons for the different EOSs. The data are as in Fig. |^. 

Fig. |[ Estimate of the background to the dielectron spectrum from the meson decays 
after freeze-out as calculated from the hydrodynamical result for the EOS A. The kinematic 



cuts and the detector resolution of the CERES experiment |jT^ are incorporated. 



Fig. 10. Dielectron spectra compared with the measurement of the CERES collaboration 



14 1 . Kinematic cuts and detector resolution are incorporated, (a) Contributions during the 
lifetime of the fireball, (b) Total dielectron spectrum including our background estimate for 
the different equations of state. 
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